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Abstract. 

We consider the Potts model in a magnetic field on an arbitrary graph G. Using a 
formula of F. Y. Wu for the partition function Z of this model as a sum over spanning 
subgraphs of G, we prove some properties of Z concerning factorization, monotonicity, 
and zeros. A generalization of the Tutte polynomial is presented that corresponds 
to this partition function. In this context we formulate and discuss two weighted 
graph-coloring problems. We also give a general structural result for Z for cyclic strip 
graphs. 
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The g-state Potts model has served as a valuable system for the study of 
phase transitions and critical phenomena pQ-[3] and has interesting connections with 
mathematical graph theory [I]-[6]. On a lattice, or, more generally, on a graph G, at 
temperature T and in an external magnetic field H, this model is defined by the partition 
function 

z = X>->" (1) 

with the Hamiltonian 

n y2 ■ i ' J (S 'K o, - u X/".- 1 > ( 2 ) 

(ij) i 

where cr, = 1, g are classical spin variables on each vertex (site) i £ G; (3 — (A^T) -1 , 
(ij) denote pairs of adjacent vertices, and Jjj are the associated spin-spin couplings. 
The graph G = G(V, E) is defined by its vertex set V and its edge (bond) set E; we 
denote the number of vertices of G as n = n(G) and the number of edges of G as e(G). 
With no loss of generality, we take G to be connected and the external field to pick 
out the value <jj = 1 from the q possible values. We first consider the case of a single 
spin-spin coupling Jjj = J and use the notation 

K = f3 J , h = PH, y = e K , v = y - 1 , w = e h . (3) 

From ([7]) below, it follows that Z is a polynomial in q, v, and w, so we write 
Z = Z(G, q, v, w) and, for the zero-field {w = 1) case, we set Z(G, q, v) = Z(G, q, v, 1). 
Positive H gives a weighting that favors spin configurations in which a^s have the value 
1, while negative H disfavors such configurations. In the limit h — > — oo, configurations 
in which any dj = 1 make no contribution to Z, so that the model reduces to the 
zero- field case with q replaced by q — 1: 

Z(G,q,v,0) = Z(G,q-l,v,l) . (4) 

The original definition of the Potts model, ([1]) and ([2]), requires q to be a positive 
integer, q G N+. This restriction is removed for the zero- field Potts model by the 
Fortuin-Kasteleyn representation [7j 

Z(G,q,v)= v< G 'k k{G,) , (5) 
G'CG 

where G' = (V, E'), E' C E is a spanning subgraph of G, and k(G') denotes the number 
of (connected) components of G'. Eq. ([5]) has the crucial property that Z{G,q,v) is 
expressed in a manner that does not make any explicit reference to the spins {crj} or 
summation over spin configurations, ■}■ This enables one to define the zero-field 
model Potts model partition function for any real q > 0. For the ferromagnetic case, 
v > 0, so Z(G, q, v ) > for q > and hence ([5]) defines a Gibbs measure. For the 
antiferromagnetic case, since v is negative (— 1 < v < 0), (J5]) does not, in general, 
yield a positive-definite Z with Gibbs measure if q ^ N + . Eq. (J5j) also establishes the 
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equivalence of the zero-field Potts partition function to the Tutte polynomial T(G, x, y), 
a function of major importance in graph theory, 

T(G,x,y)= £ (* - l)*<^-*<% - 1)*^ , (6) 

G'CG 

where c(G') = e(G') + k(G')—n(G') is the number of independent cycles on G' pf]-[6],[E]- 
[lig. The equivalence is Z(G,q,v) = (x-l) k ^{y-l) n ^T{G, x, y), where x = l+(q/v). 

The Fortuin-Kasteleyn cluster formula (jSJ) was generalized to the case of a nonzero 
external magnetic field by F. Y. Wu [2| [TT]. Denote each of the connected components 
of G' as G'i, i = 1, fc(G'). Wu's result is [21 [UJ 

KG') 

Z(G, q, v,w)=Y, ^ II (9 " 1 + ^ n(G °) ' ( 7 ) 

G'CG i=l 

We first use the Wu formula d^J) to prove a number of properties of g, t>, to) 
concerning factorization, monotonicity, and zeros. Combining (jHJ), which shows that 
Z(G,q,v) contains a factor of q, with we deduce that Z(G, q, v,0) contains a 
factor of (g — 1). Substituting g = in ([7|) and using the factorization io n< - G ^ — 1 = 
to Y^l=o + where to = to — 1, we prove that Z(G, 0, t>, 10) contains a factor of 

(to — 1). Setting to = g — 1 in (ED yields the result that Z(G, q,v,q — 1) has (g — 1) as a 
factor. Substituting w = in ([7j) is another way to derive (j4j). Two elementary results 
are Z(G, 1, v, w) = (v + \)<G) w n{G) and Z (q^ g? q, = (g - 1 + w)< G \ 
We can write Z(G,q,v,w) in several equivalent ways: 

n(G) e(G) n(G) e(G) 

z(G,q,v,w) = ^2 a rstq r v s w t = ^2^2 b rs tq r y s w t 

r,t=0 s=0 r,t=0 s=0 

n(G) e(G) n(G) e(G) 

= ^ Yl c rstq r v s w t = drstq^'w 1 , (8) 

r,t=0 s=0 r,t=0 s=0 

where w = w — 1 as before, g = g — 1, and the coefficients a rst , b rst , c rst , and d rst are 
integers. Some a rst and 6 rsi can be negative. In contrast, the Wu formula (jTj) shows 
that all of the nonzero c rs t are positive. This leads to three monotonicity and zero-free 
properties in the corresponding variables q, v, and w, taken here as real: (i) for q > 
and v > 0, Z(G, q, v, w) = Z is a monotonically increasing function (MIF) of to > and 
has no zeros on the positive w axis; (ii) for v > and to > 0, Z is a MIF of g > and has 
no zeros on the positive g axis; (iii) for w > and g > 0, Z is a MIF of v > and has no 
zeros on the positive v axis. We can also prove that all of the nonzero d rst are positive 
by using (JTj) together with the relation used above, w 11 ^ — 1 = to Y^=q 1 (^ + 
Since each term in the expansion of (to + 1) is positive for each G[, this shows that 
the nonzero d TSt are positive. This yields three more monotonicity and zero-free results 
(which have some overlap with (i)-(iii)): (iv) for q > and v > 0, Z is a MIF of to > 
and has no zeros on the positive to axis; (v) for v > and to > 0, Z is a MIF of g > 
and has no zeros on the positive g axis; (vi) for to > and g > 0, Z is a MIF of f > 
and has no zeros on the positive v axis. Monotonicity relations for borderline cases are 
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covered by our results above; e.g. for q = 1, Z(G, 1, v, w) is a MIF of v > — 1 for w > 
and a MIF of w > for v > — 1; for u = 0, g, 0, w) is a MIF of q - 1 + w > 0, etc. 
We define a rational function that generalizes the Tutte polynomial, namely 

U(G,x,y,w) = (x-l)- k(G \y-iy n ^J2(y- 1 ) e{G ' ) x 

G'CG 

KG 1 ) 

x Yl ( x v - x - y + w n{G,i) ) ■ (9) 

8=1 

This function satisfies U(G,x,y,w) = (x — l)~ k ( G '(y — l)^ n ^ Z(G, q, v , w) and 
U(G, x,y, 1) = T(G,x,y). Although T(G,x,y) and Z(G,q,v) satisfy deletion- 
contraction relations, we note that for w not equal to 1 or 0, the functions U(G, x, y, w) 
and Z(G,q,v,w) do not, in general, satisfy such deletion-contraction relations. 

We define two types of graph coloring problems and use special cases of (J7|) to 
describe these. Although graph coloring has been investigated intensively [l]-[6], [9j- 
[TU] . [12] . these two types of graph colorings have not, to our knowledge, been studied 
before. Recall that the chromatic polynomial P(G, q) counts the number of ways of 
assigning q colors to the vertices of a graph G such that no adjacent vertices have the 
same color. This "proper g-coloring" of the vertices of G is equivalent to Z for the 
zero-temperature, zero-field Potts antiferromagnet, v = —1: P(G,q) = Z(G,q, —1). We 
generalize this to a weighted proper g-coloring of the vertices of G, as described by 
the polynomial Ph(G,q,w) = Z(G, q, — 1, w). For H < 0, i.e., < w < 1, we have 
a weighted graph coloring problem in which one carries out a proper g-coloring of the 
vertices of G but with a penalty factor of w for each vertex assigned the color 1. For 
H > 0, we have a second type of weighted graph coloring problem, namely a proper 
vertex g-coloring with a weighting that favors one color. Since this favoring of one color 
conflicts with the constraint that no two adjacent vertices have the same color, the range 
w > 1 involves competing interactions and frustration. 

Both of these weighted graph coloring problems have physical applications. For 
example, the weighted coloring problem with < w < 1 describes the assignment 
of frequencies to commercial radio broadcasting stations in an area such that (i) 
adjacent stations must use different frequencies to avoid interference and (ii) stations 
prefer to avoid transmitting on one particular frequency, e.g., because it is used for 
data-taking by a nearby radio astronomy antenna. The graph coloring problem with 
w > 1 describes this frequency assignment process with a preference for one of the q 
frequencies, e.g., because it is most free of interference. We note some other special 
cases. Just as the Tutte polynomial T(G, 1 — q, 0) gives, up to a prefactor, P(G, q), so 
also T(G, 0, 1 — q) determines the flow polynomial F(G,q), which counts the number 
of nowhere-zero g-flows on G that satisfy flow conservation mod q at each vertex. The 
function U(G,0, 1 — q, w) then defines a weighted flow problem. With < w < 1, this 
could describe a discretized flow analysis in an electrical circuit or traffic flow situation 
in which one incorporates a finite penalization for one, say the maximal, flow, in order 
to minimize power dissipation in resistors in the circuit case or to minimize traffic jams 
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^ = E(- 1 ) i ( 2d j %- l ) d ' j - en) 



in the traffic case. 

For a planar G, P(G, q) counts not just the number of proper g-colorings of G 
vertices but also, equivalently, the number of proper g-colorings of the faces of the dual 
graph G*. Similarly, for planar G, Ph(G,q,w) is a measure not only of the weighted 
proper g-colorings of G vertices, but also, equivalently, the weighted proper g-colorings 
G* faces. 

We have used (j7]) and combinatoric arguments of the type in (13] to obtain a 
general structural determination of Z(G,q,v,w) for cyclic and Mobius strip graphs G s 
of regular lattices of fixed width L y vertices and arbitrary length as well as self-dual 
strips of the square lattice, extending [14J. This length is L x = m (L x = 2m) for square 
and triangular (honeycomb) strips. For cyclic G s we find 

L y n Z h(L y ,d) 

Z(G s ,q,v,w) = J2 E ^ d \q)[X Gs , LyAj (q,v,w)] m , (10) 
where Zh connotes Z for h ^ and 

d 
3=0 

We have nzh(L y , L y ) = 1, 72^^(1,0) = 2, and nzh(L yi d) = for d > L y \ the other 
nzh(L y , d) are determined by the recursion relations n Z h(L y + 1,0) = 2n Z h(L y , 0) + 
nzh{L y , 1) and, for 1 < d < L y + 1, 

n Z h{L y + 1, d) = n Z h{L y , d-l) + 3n Z h{L y , d) + n Z h{L y , d + 1) . (12) 

The form for Mobius strips involves switches of certain c^'s (generalizing switchings 
in the w = 1 case [13]), which are given in detail elsewhere [15]. For these cyclic 
(and Mobius) strip graphs of width L y , the total number of different A's, N Z h,L y = 
J2dlo n zh(L y ,d), is 

It is straightforward to generalize ([7]) to the case where the spin-spin couplings J^- 
depend on the edges e#. Let us define Kij = j3Jij, yij = e Kij , Vij = yij — 1, and the set 
of v e for e = eij G E as {v e }. Then we have 

KG') 

z(g, g , {v e }, w) = ^2 [ n v ] [ n (« - 1 + wn(G '^] ■ ( 14 ) 

G'CG eeE' i=l 
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